%% 5.2	Soil physics

close all; clearvars; clc;

%% 5.2.a
% Calculate the porosity of a dry soil if it is known that the bulk density
% is 1.8 [g/cm3] and that the soil exists of sand and silt.

% Assuming all grains are made from quartz
rho_quartz  = 2650;     % [kg/m3]
rho_bulk    = 1800;     % [kg/m3]

fprintf('The porosity can be estimated at %f [-]\n',1-rho_bulk/rho_quartz);


%% 5.2.b
% Calculate the dry soil bulk density if the soil consists of fine sand 
% with a porosity of 35%.

p   = 0.35;     % [-]

fprintf('The density of the bulk soil is %f [kg/m3]\n',(1-p)*rho_quartz);
